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Abstract. In this paper we investigate overdetermined systems 
of scalar PDEs on the plane with one common characteristic, whose 
general solution depends on 1 function of 1 variable. We describe 
linearization of such systems and their integration via Laplace 
transformation, relating this to Lie's integration theorem and for- 
mal theory of PDEs. 



Introduction 

Consider an overdetermined system £ of partial differential equa- 
tions, which we assume formally integrable (i.e. all the compatibility 
conditions fulfil) and regular (this is a generic condition micro- locally). 
We will restrict to systems with two independent and one dependent 
variables (but the theory can be extended to other systems). 

The characteristic variety Char c (£) of this system is an effective 
divisor on CP 1 , i.e. is a collection of points with positive multiplicities. 
Let uo = deg(Char {£)) be the total multiplicity. This number was 
called class of the system £ by Sophus Lie |L]] . 

In terminology of Ellie Cartan uo is the Cartan integer s\ (provided 
the Cartan character is 1: S2 = 0). Cartan's test |CJ 2 | implies that the 
general solution u G Sol(£) depends on uo functions of 1 variable. 

Note that 00 can be described in a different way: Since the character- 
istic variety is discrete, the symbol gt of £ stabilizes: lim dim^ = 00 

k—^oo 

(when the system becomes involutive, see Appendix [B] for details). 

The case uo = corresponds to finite type systems and integration 
of £ can be reduced to a system of ODEs via the Frobenius theorem. 

Another well-known class is u; = 2, especially scalar second order 
PDEs on the plane. One of the classical approaches to such systems is 
the Laplace transformation. 

In this paper we show that Laplace transformations exist in the case 
uo = 1 as well (a comment on the case of general uo will be made at the 
end of the paper). 



Key words and phrases. Lie's class 1, overdetermined systems of PDEs, charac- 
teristic, integral, symbols, compatibility, Spencer cohomology. 
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0.1. Main results. In his paper of 1895 [Li] S. Lie demonstrated that 
a compatible (=formally integrable) class u — 1 overdetermined system 
is integrable by reduction to ODEs. Modern proof and applications of 
this result will be discussed in |K 2 | . 

In this paper we demonstrate that in certain cases this reduction 
is very precise, and can be decomposed into a sequence of external 
transformations in jets that are analogs of the classical Laplace trans- 
formation. These latter are the differential substitutions with differ- 
ential inverses. This provides a method for an effective (algorithmic) 
integration of PDE systems of class uj = 1. 

The transformations will be fully described in the linear case, and 
linearizable systems will be characterized via a simple criterion. For 
non-linear systems we discuss some phenomenology and examples. 

Assume that the only characteristic for u = 1 linear system is 
straightened X = d x (this involves integration of a 1st order non- 
autonomous non-linear scalar ODE). 

Theorem 1. A linear system £ of class uj = 1 with one characteristic 
X is integrable in closed form and quadratures for generic £, or is 
reducible to class uj = (ODEs) in singular cases. 

Thus a generic linear system of class uj = 1 has representation for 
the general solution via a differential operator and mild nonlocality 
(quadrature) applied to an arbitrary function. This is neither true for 
for uj > 2 (e.g. scalar 2nd order PDEs on the plane), nor for uj = 
(e.g. linear 2nd order ODEs, equivalent to Riccati equation). 

Remark. This distinguishes uj = 1 linear systems, but does not extend 



to nonlinear uj = 1 class systems, which will be discussed in K 2 



Internal geometry of linear/linearizable systems is quite simple: they 
correspond to Goursat distributions with growth vector (2,3,4,. . . ). On 
the other hand the external geometry (which is governed by the pseu- 
dogroup of point triangular transformations) is rich and is character- 
ized by differential invariants ( |Ci| contains an example of two sys- 
tems which are equivalent internally but not externally, they have type 
E 2 + E 3 in the notations of §|2J E.Cartan proved that this is not possible 
for 2E2). Some of these appear naturally in our approach, and are the 
analogs of the classical Laplace invariants. 

In this paper we characterize class uj = 1 systems from the external 
viewpoint, and the most important invariant will be the complexity x 
of the system £ . This is an integer- valued non-negative quantity which 
always decreases under Laplace transformation. Some other quantities 
characterizing the size of the solutions space Sol(£ ) are also relevant. 
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In the future work |K 2 | , we shall model reduction types, based on 
the normal forms of rank 2 distributions. Cartan-Hilbert equation is a 
classical example of this general construction for class u = 1. 



0.2. Background and outline. We will exploit the geometric theory 
of PDEs, namely jet-geometry and Spencer formal theory [S]. The 
reader is invited to consult |KLj| or a short exposition in Appendix |B] 

Also the geometry of distributions will be occasionally used in a 
minor part of the text. The structure of the paper is as follows. 

In Section [TJ we present a study of class u = 1 linear systems with 
low-complexity. We introduce the pseudogroup of frame/coordinate 
changes and calculate some relative invariants. Then we define gener- 
alized Laplace transformations on phenomenological level as differential 
substitutions, which simplify the system (and have differential inverses 
modulo the equation). 

This latter condition means that the Laplace transformations de- 
crease the complexity — the notion that is introduced in Section [2j 
There we also describe the totality (zoo) of all systems of class u = 1 
and discuss properties of the complexity. 

The generalized Laplace transformations are rigorously defined in 
Sectional Then we prove the main results about existence, uniqueness 
and effectiveness of Laplace transformations for the linear case. 

Section H] discusses some features of the non-linear situation. Here 
new ideas are required: integrable extensions, non-Moutard form of 
solutions etc. Some of these will be discussed in forthcoming paper 



K2I , so we restrict to several examples. 

In Section we present a short historical overview and give a brief 
discussion of possibilities and difficulties of generalizations of the theory 
for systems of class u> > 1. 

The appendices supplement some of the background material. 



Acknowledgment. It is a pleasure to thank Nail Ibragimov for his 
translation of Sophus Lie paper [L]] in |LG] , which was a starting point 
for this paper. I am grateful to Valentin Lychagin for many discussions 
on the initial stage of the project. The paper is strongly influenced by a 
collaboration with Ian Anderson, to whom I am thankful. Hospitality 
of MLI (Stockholm) in 2007, MFO (Oberwolfach) in 2008, Banach 
center (Warsaw) and Utah State University in 2009, as well as IHES 
in 2010 is acknowledged. 



4 



BORIS KRUGLIKOV 



1. Linear systems on plane with one simple characteristic 

Consider a compatible PDE system £ of class u — 1, meaning there 
is only one characteristic^] of multiplicity 1 (it can be taken X = d x for 
linear and some quasi-linear systems, but not for a fully non-linear £). 

The simplest situation is one 1st order PDE, which is classically 
known to be solvable via the method of characteristics. In general £ is 
generated by r PDEs of possibly different orders and we list only gen- 
erators (disregard prolongations) of the system. The number of these 
generators can be determined invariantly via the Spencer 5-cohomology 
groups: r = dim if*' 1 (£). 

Starting from some jet-level we have for the symbol dimc^ = 1. This 
imposes restriction on the form of the system and we shall investigate 
them successively according to complicacy, which includes the number 
r, orders of the PDEs and orders of the compatibility conditions. 

Spencer cohomology have dimensions h° = dim if*' = 1, h 1 = 
dim if*' 1 = r and from vanishing of the Euler characteristic]! we find 
h 2 = dim if*' 2 = h 1 — h° = r — 1. This is the number of compatibility 
conditions our system £ satisfies. 

In this section we consider only linear PDEs in one (scalar) unknown 
u = u(x, y), in which case the characteristic is a vector field on the base 
M = M. 2 . We will denote it by X and its complement by Y. 

We take the freedom of considering X, Y to be first order differential 
operators, rather than vector fields (and still call it frame), and the 
symbol of X is the characteristic. 

The maximal transformation group Q of this nonholonomic frame is 

X ^ KX + a, Y h-> qY + bX + c. 

We can always arrange [X, Y] = and in doing sell we get functional 
dimension of the pseudogroup over the algebra A = C°°(M) equal 2. 

We however relax the requirement to ordLY, Y] = (this is always 
achieved by a choice of n, q) and the pseudogroup with this condition 
has functional dimension func.dim^C?) = 3. 

Generators of the pseudogroup Q are the following transformations 

I^I + a, y^F (basic gauge) 

X i — y e x X, Y \-t Y + bX + c (auxiliary gauge) 

X i-> kX, Y i — y qY (auxiliary gauge) 

^he characteristic covectors are defined as in the formal theory of PDEs [Sj 
|KLi| , see Appendix [Bl For 2-dimensional base they dualize to characteristic vectors 
(defined up to non-zero multiple). 

2 This holds for all systems, reducible to ODEs (in our case oj = 0, 1). 

3 With this the pseudogroup is equivalent to the linear-triangular pseudogroup 
x i — y x(x, y), y i — y 7(2/), u 1— > -&(x, y)u provided the characteristic is d x . 
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with a,b,c e C7°°(M) arbitrary and X,K,q e C°°(M) satisfy [Y + 
bX, A] = [X, b], [X,q] = \Y,k) = 0, i.e. A, K, ? are functions of one 
argument (we write [X, 6] instead of X(b) as X is a 1st order differential 
operator and similarly in other cases). 

Remark. Using bigger pseudogroup Q instead of point transformations 
is harder from the viewpoint of differential invariants, but is more con- 
venient from factorization viewpoint that is our main goal now. 

Note that alternatively we could use the pseudogroup Q' : X i— > X + a, 
Y i— ► Y + bX + c (which satisfies the condition [X, Y] — mod X + 
order 0) with func.dim^t?') = 3, achieving the same results. 

Notice that the basic gauge transformations form a normal subgroup 

H, the auxiliary form a subgroup /C and we get H ■ /C = Q. We'll be 
using the following easy result: 

Proposition 2. Let Q — 1~L ■ K be a decomposition of a group into 
the product of two subgroups the first of which is normal. Consider an 
action of Q on a space M (manifold in finite- dimensional case) and 
suppose there exists a global H-transversal subspace L C M invariant 
with respect to JC. Then the space of invariant functions C°°{L) pulled 
back to M via %-action coincides with the space of Q -invariants: 

C°°(M) e ~ C°°{L) K . 

Now we consider low complexity linear systems of u — 1 type to 
illustrate the general picture that we'll sum up later. 

I. 1. Two equations of the 2nd order: 2E2. Consider the case of 
two PDEs (r = 2) on the plane M. As the common characteristic ia^ 
X, the equations have symbols X 2 and YX (it will be more convenient 
to write this order for calculations later on). Thus our system S is 

J X 2 u + a x Xu + b{Yu + C\U = 0, 
\ YXu + a 2 Xu + b 2 Yu + c 2 u = 0. 

This system has one compatibility condition since h 2 = 1. It implies 
bi = 0. By the basic gauge (absorbing functions into X) we transform 
the system to have b 2 = 0. Then compatibility implies c\ = 0. 

Namely we have X 1— >■ X = X + b 2 , and for the transformed 2E 2 we 
have b 2 = c\ = and 

c 2 = c 2 - a 2 b 2 - [Y, b 2 ] 

(compare with the classical Laplace invariants, see Appendix |A|). This 
c 2 is a relative differential invariant with respect to auxiliary gauge, 
and so is a relative invariant of the pseudogroup Q action (here and in 
what follows we use Proposition [2]). 



We write just X instead of more appropriate symb(AT). 
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Two sub-cases are possible (we are omitting tildes from now on). 

Y| 2 : c 2 7^ 0. Then differential substitution^ v = Xu leads us to one 
first order equation (essentially parametrized ODE): 

(X + a x )v = 0. (1) 

The second PDE in £ yields the inversion formula u = — c 2 " 1 (F + a 2 )v. 
Note that in this case the inversion is differential. 

The equation [Xc 2 l (Y + a 2 ) + l]v = comes from composition of 
the substitution and the inversion, but it is a differential corollary of 
dlj) due to compatibility of £. 

Y 22 : c 2 = 0. Then differential substitution v = Xu leads us to the 
Frobenius system: 

(X + a x )v = 0, (Y + a 2 )v = 0, (2) 

which is compatible (due to compatibility of £) and has 1-dimensional 
vector space of smooth solutions. 

However inversion in this case is integral u = X _1 v (note however 
that this is a parametrized ODE) providing one unknown function for 
the general solution u via essentially one solution of (|2J). 

1.2. Two equations of the 2nd and 3rd orders: E 2 + E 3 . In this 
case (r = 2) a common characteristic condition leads us to two cases. 
Since we consider first the lower order PDEs, prolong them, add new 
equations and so on, the general position case will be such that the 
lower order symbols are in general position. The first of the cases 
below is of general position, the second as we shall see has degeneracy. 

T 23 : The symbols of the PDEs are YX and X 3 . Thus the system 
£ has the form: 

f YXu + dXu + diYu + e x u = 0, 

\ X 3 u + a 2 X 2 u + b 2 Y 2 u + c 2 Xu + d 2 Yu + e 2 u = 0. 

The compatibility condition yields b 2 = d 2 = and with the basic 
gauge we get d x = 0. Then compatibility leads to e 2 = 0. 
After this transformation the new coefficient e x (which is 

e?i = ei — c\d\ — [Y, d x ] 

via the old coefficients and we remove tilde; such precise formulae will 
be omitted in the future cases) is a relative invariant with respect to 
auxiliary gauges (here the pseudogroup Q is smaller as the symbol of 
£ fixes directions of both X and Y thus reducing some of auxiliary 
gauges), and consequently with respect to Q. 
The following sub-cases occur: 



'Such substitutions, simplifying £ , will be called Laplace transformations. 
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Y23 : ei 7^ 0. The substitution Xu = v has differential inverse from 
the first equation: u = —e±(Y + C\)v. 

The second PDE of £ leads to a 2nd order PDE for v with symbol 
X 2 . Inserting the inversion into the substitution we get another PDE 
with symbol YX. The obtained system 2E 2 is: 

f (X 2 + a 2 X + c 2 )v = 0, 
\ (Xe^(Y + Cl ) + l)v = 0. 

Y23 : e\ = 0. Then the substitution t> = X-u yields two PDEs, one of 
the 1st and the other of the 2nd order, but this system is of Frobenius 
type (and compatible): 

(Y + Cl )v = 0, (X 2 + a 2 X + c 2 )v = 0. 

The inversion is however integral: u = X~ 1 v. 

Y23 : The symbols of the PDEs are X 2 and Y 2 X. Thus the system 
£ has the form: 

J X 2 u + c\Xu + d{Yu + e\u = 0, 

\ Y 2 Xu + a 2 YXu + b 2 Y 2 u + c 2 Xu + d 2 Yu + e 2 u = 0. 

Compatibility condition yields d\ — and with the basic gauge we 
obtain 62 = 0. Then the compatibility leads to e% = 0. 

The new coefficient d 2 is a relative differential invariant and e 2 is a 
relative invariant provided d 2 = 0. The following are the sub-cases: 

Y23 : d 2 7^ 0. The substitution v = Xu to the first PDE yields the 
first order equation 

(X + Cl )v = 0. (3) 

But the inversion is given by the solution of the compatible Frobenius 
system formed by the substitution and the second PDE: 

Xu =v, (Y + d 2 l e 2 )u = -d^iY 2 + a 2 Y + c 2 )v. 

Let us remark that integration of this system adds not only one con- 
stant, but also one derivative to the arbitrary function of 1 variable 
coming from the general solution of ()3]) . 

Y23 : d 2 = 0, e 2 7^ 0. Here v = Xu has the differential inversion 
u = —e 2 1 (Y 2 +a 2 Y+c 2 )v. First order PDE ([3]) remains the transformed 
equation with this substitution. 

It seems that insertion of the inversion to the substitution leads to a 
3rd order PDE with the symbol Y 2 X, but it is a differential corollary 
of 03]) due to compatibility of £. Thus we get only one first order PDE. 

Y23 : d 2 = 0, e 2 = 0. Here substitution v = Xu leads to the Frobe- 
nius system of the 1st and 2nd order equations: 

(X + Cl )v = 0, (Y 2 + a 2 Y + c 2 )v = 0. 
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The inversion operator u = X v is integral. 

1.3. Three equations of the 3rd order: 3E3. In this case (r = 3) 
S has one single characteristic X iff its form is: 

X 3 u + a x X 2 u + b x YXu + c x Y 2 u + d x Xu + e x Yu + f x u = 0, 
YX 2 u + a 2 X 2 u + b 2 YXu + c 2 Y 2 u + d 2 Xu + e 2 Yu + f 2 u = 
F 2 Xm + a 3 X 2 u + 6 3 FXu + c 3 Y 2 u + d 3 Xu + e 3 Yu + f 3 u = 0. 

Compatibility implies c x = c 2 = and the basic gauge yields c 3 = 0. 
With this the compatibility gives: e x = e 2 = 0. 

The new coefficients e 3 , b±, f\ are relative differential invariants (un- 
der auxiliary gauge and so under Q), as well as /^(upon f\ = 0), 
/ 3 (upon fi = f 2 = 0) (which are conditional invariants). 

Compatibility ties them so: f\ = &ie 3 , f 2 = b 2 e 3 + [X, e 3 ]. The 
following are the cases: 

T^ 33 : e 3 ^ 0,6i ^ & fi ± 0. The fact that 61 ^ implies 
that the first two equations are not compatibly in itself, and that 
their compatibility condition implies the third PDE of £. Thus we can 
discard the latter. 

Substitution v = Xu can be inverted from the first PDE: u = 
-fr\X 2 + a x X + b x Y + d x )v. 

A linear combination of the second and the first PDEs from 8 yields 
a 2nd order equation on v. Another equation is of the third order and 
is obtained by inserting the inversion into the substitution. Thus we 
get a system on v of type E 2 + E 3 (T^): 

f ((YX + a 2 X + b 2 Y + d 2 ) - f x l f 2 (X 2 + a x X + b x Y + d x )) v = 0, 
\ (Xf^iX 2 + a x X + b x Y + d x ) + l)v = 0. 

The symbols of these equations after an auxiliary gauge are YX, X 3 
and the system is compatible. 

T333 : e 3 7^ 0, f x = =3- b x = 0. Suppose first f 2 ^ 0. 

Substitution v = Xu can be inverted from the second PDE: u = 
-L[v], where L = f 2 \YX + a 2 X + b 2 Y + d 2 ). 

Then the first PDE of £ and insertion of the inversion to the third 
PDE yield the compatible system on v of type E 2 + E 3 (T| 3 ): 

j (X 2 + a x X + b x Y + d x )v = 0, 

1 ((e 3 y + h)L - {Y 2 + a 3 X + b 3 Y + d 3 ))v = 0. 

T333 : e 3 7^ 0, f x = 0, b x = 0, f 2 = 0. The substitution v = Xu turns 
the first two PDEs of 8 into 2E 2 : 

(X 2 + a x X + d x )v = 0, {YX + a 2 X + b 2 Y + d 2 )v = 0, 

6 The term " compatible" means " formally integrable" |KLi| , and is more restric- 
tive than (formally) " solvable" , see Appendix [Bj 
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with the inversion found from the compatible system of Frobenius type: 
Xu = v, (Y + e^f 3 )u = -e^(Y 2 + a 3 X + b 3 Y + d 3 )v. 

Y333 : e 3 = 0, implying fi = f2 = 0. Further exploration of compat- 
ibility implies 61/3 = 0. Assume first f 3 ^ =r- b\ = 0. 

Then the third PDE of E implies inverse u = — f 3 ~ 1 (Y 2 + a 3 X + b 3 Y + 
d 3 )v for the substitution v = Xu. The first two PDEs of S yield the 
same 2E 2 as in T333. 

f 333 : e 3 = 0) /1 = h = fs = 0- Then substitution v = Xu leads to 
the Frobenius system 



(YX + a 2 X + b 2 Y + d 2 )v = 0, (Y 2 + a 3 X + b 3 Y + d 3 )v = 0. 



with the integral inverse u = X 1 v. 

1.4. Two equations of the 3rd order: 2E 3 . This case (r = 2) is 

a bit more complicated since the compatibility conditions are given by 
operators of the 2nd order. The characteristic will again be denoted 
by X. Thus the two symbols have the form X ■ Qi,X ■ Q 2 , where Qi 
are linearly independent quadrics on T*M, i.e. elements of S 2 TM. 

Denote by IT 2 = (Q U Q 2 ) the plane in S 2 TM. Let q e S 2 T*M be 
a non-zero element of Ann(II 2 ) and coy € T*M a non-zero element of 
Ann(X). The quadric q on TM is not a square, since otherwise Q\ and 
Q 2 would have a common factor, which would be another characteristic. 

Depending on the position of q regarding the cone of degenerate 
quadrics and the line u Y we obtain the following two different possi- 
bilities (classification over C is by one case smaller) with ux being a 
covector, complementary to uy- 



The corresponding normal forms for the plane n 2 generate the normal 
forms of the symbol of S. They fix the direction of Y, and the auxiliary 
group becomes smaller (reduction of Q similar to the case T23). 

Y33 : The symbols of the two equations are (X 2 ± Y 2 )X, YX 2 . The 
equation has the form 



YX 2 u + a x X 2 u + b x YXu + c x Y 2 u + d x Xu + e{Yu + f x u = 0, 



(X 2 ± Y 2 )Xu + a 2 X 2 u + b 2 YXu + c 2 Y 2 u + d 2 Xu + e 2 Yu + f 2 u = 0. 



Compatibility implies that c\ = 0. The basic gauge yields c 2 = and 
then compatibility gives e± = 0. 

Then the new functions e 2 , &i, f\ are relative invariants and the com- 
patibility ties them as f\ = b\t 2 \ f 2 is a conditional invariant on e 2 = 0. 




q = co x =F ijjy or q — w^uy. 
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T33 : fx 7^ =>- e 2 7^ 0. Substituting Xu = v we get the differential 
inversion u = —Li[v], L\ = (YX + aiX + b{Y + di), which be- 
ing inserted into the second PDE of £ and the substitution gives two 
differential equations, with the notation L 2 = X 2 ±Y 2 + a 2 X + b 2 Y + d 2 : 

f ((e 2 Y + f 1 )L 1 + L 2 )v = 0, 

\ (XL 1 + l)v = 0. W 

This system is however not compatible. Indeed, we obtained it with 
the help of inversion from the system Xu = v, (e 2 Y + f 2 )u = —L 2 [v}. 

Writing the compatibility condition for this system, i.e. multiplying 
from the left by (Y + pi) and by (—X + p 2 ) respectively and adding, 
we obtain an equation of the type 9 ■ u = L 3 [v], where L 3 = XL 2 + . . . 
and 9 7^ because the compatibility conditions for £ are given via 
differential operators of the second (not by first!) order. 

Combining the latter equation with the PDE u = —Li[v] we get 
a 3rd order differential equation on v with the symbol (X 2 ± Y 2 )X, 
which is the compatibility condition for system (j3J). Finally taking 
linear combinations of the derived PDEs we obtain the system of type 
3i?3 on v with the symbols X 3 , YX 2 , Y 2 X, which is now compatible 
due to compatibility of £. 

T33 : fi = 0, e 2 7^ 0. Then substitution Xu = v together with (e 2 Y+ 
f 2 )u = —L 2 [v] has Frobenius type and its compatibility condition gives 
a 3rd order PDE on v. Another PDE on v is given by the first PDE of 
£ and it is has 2nd order. This pair E 2 + E 3 on v has symbols YX, X 3 
(T23) and is compatible. The inversion is however Frobenius. 

T33 : j\ = 0, e 2 = 0, f 2 7^ 0. Then we can get differential inversion 
from the second PDE of £ u = —f 2 L 2 [v]. The first PDE of £ becomes 
the 2nd order PDE (YX + a 2 X + b 2 Y + d 2 )v = 0, while substitution of 
the inversion into Xu = v yields (Xf 2 l L 2 + l)(f) = 0, which modulo 
the previous PDE is a 3rd order equation with symbol X 3 . This is a 
compatible system of type E 2 + E 3 (T| 3 ). 

T33 : fi — 0, e 2 = f 2 — 0. The substitution v = Xu reduces the sys- 
tem to a pair of compatible 2nd order PDEs with symbols X 2 ±Y 2 , YX. 
Hence the system is of finite (Frobenius) type and the substitution has 
integral inverse. 

T33 : The symbols of the two equations are X 3 ,Y 2 X and so the 
system has the form 

f X 3 u + a x X 2 u + biYXu + c x Y 2 u + d x Xu + e x Yu + f\u = 0, 
\ Y 2 Xu + a 2 X 2 u + b 2 YXu + c 2 Y 2 u + d 2 Xu + e 2 Yu + f 2 u = 0. 

By compatibility c\ = 0, basic gauge gives c 2 = 0, then compatibility 
implies e\ — 0. 
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Then the new functions e 2 , b±, f\ are relative invariants and the com- 
patibility ties them as f\ = b\e 2 \ f 2 is a conditional invariant on e 2 = 0. 
Let Li = X 2 + ai X + bxY + d u L 2 = Y 2 + a 2 X + b 2 Y + d 2 . 

T33 : fi 7^ 0. The inversion is differential u = — f± Li[v] (from 1st 
PDE of £ ). Insertion of this into the 2nd PDE of £ and into substitution 
v = Xu, together with addition of compatibility conditions to the last 
two leads to a compatible system 3E3 with symbols X 3 , YX 2 , Y 2 X. 

"""33 : /1 = 0) e 2 7^ 0. The first PDE on v is L\[v\ = 0. The inversion 
can be found from the system 

Xu = v, (Y + e 2 l f 2 )u = -e 2 l L 2 v. (5) 

The compatibility of this system has the form [Xe 2 x L 2 + . . . )v = ku. 

In the case k = we get the 3rd order PDE on v with symbol Y 2 X, 
so the reduced system is E 2 + E 3 (T| 3 ). The inversion is a Frobenius 
system. 

If k ^ 0, then the inversion is differential u = k~ l (Xe 2 L 2 + . . . )v . 
Substituting this into the 2nd PDE of (JSJ) yields an equation of 4th 
order with symbol Y 3 X. Thus the reduction is of type E 2 + E±. 

*"*33 : /1 = e 2 = 0, f 2 ytz 0. The inversion is differential u = 
— f 2 L 2 [v]. The first PDE on v is Li[v] = and the second is ob- 
tained by insertion of the inversion into substitution Xu = v. The 
symbols of the obtained equations are X 2 ,Y 2 X. They are compatible 
and have type E 2 + E 3 (T23). 

T33 : fi = 0, e 2 = 0, f 2 = 0. The reduction is a Frobenius system 
of the second order Lx[v] = 0, L 2 [v] = and the inversion is integral 
u = X~ 1 v. 

1.5. Examples of calculations. Let us illustrate integration of linear 
co = 1 systems for some representative types discussed above. 

Example 1. The following compatible system has type 3E 3 (T333). 

Uxxx — ^2 ^xx ~ t~~ ^3" ~ t~~ ^3 ^x ~£~3 ^ 

_ 4z 2 +6x+18 u ^ _± u I 8x 2 -6x+36 , 18 u 
xxy <$ X y xx x 2 xy 3x 2 y x x 2 y 

_ 2x 3 +9x 2 +45a:+54 , -5z 2 +12x+18 

x yy ~ 9p Uxx ' 3xy Ux v 

1 -16x 3 +9:r 2 -36:r-108 „, , 3 , 4x 2 -9x-18 „, 

~T~ rt 9 U*T. I |~ o LI 

\ixy £ y y %y 

It reduces by differential substitutions to the one 1st order PDE 
u x = as follows 3E3 — > E 2 + E$ — > 2E 2 — > E\. The branch is generic, 
i.e. in each type we have a generic case. 

We call the above arrows (generalized) Laplace transformations. They 
all have differential inverses. The composition of inversion formulae 



12 BORIS KRUGLIKOV 

gives the solution: 

u = 9y 3 f"\y) + 27xy 2 f"(y) + 36x 2 yf(y) + 16x 3 f(y). 



Example 2. Another compatible 3E 3 of type Y333 is 



^^xxx — ^2 ^^xx I ^2 ^~^xy rj*2. x '^'y ^ x^ 

U*XXy x 2 U XX ~ x 3 ^xy x 2 U X x g Uy ~ X Z U 

_ 4s-2y ,2x^+2xy ,3 §_ u -§y , fi 

Now the route of Laplace transformations is 3E 3 — > E 2 + E 3 — >■ 2E 2 
and it is not generic, as E 2 + E 3 has type T23, and 2E 2 has type T22, 
so this reduces to a Frobenius system of the 1st order and the inverse 
is integral. The composition of inversion formulae gives the solution: 

u = x 3 f"(y) - 6x 2 f'(y) + 6xf(y) + Cy. 



Example 3. Finally consider a system of type E 2 + E 3 , which has 
inverse Laplace transformation of the Frobenius (uu = 0) type. 

^xx ^xyy ^xy Uy. 

y y 

Laplace transformation v — u x reduces this to the equation v x = 0, but 
the inverse is given by the compatible system 

f u x = v, 

|^ Uy = XV y — yVyy. 

General solution to equation E\ on v together with integration of the 
Frobenius system yield 

u = (x + l)<p(y) - yip'(y) + C. 

Remark. In Examples 2 and 3 we displayed the "happy" cases when 
the quadratures are expressed in the closed form. Generally the solution 
of a Frobenius system cannot even be reduced to the quadratures. 



2. Classification of types ofw = 1 systems 

In this section we split the totality of u — 1 systems into classes, and 
discuss transformations between them. Most results of this section are 
general and apply equally well to non-linear uj — 1 systems. 
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2.1. Zoo for oj = 1. Let us compose a table, where we put compatible 
systems £ of class ou — 1, and organize its columns by the maximal order 
A; max and its rows by the number of equations r = h 1 = dim H*' 1 (S). 

We introduce the following rule for the choice of generators of S. 

Consider the orders of the system: /c min — k± < • • ■ < k r — k mSLX , 
and denote the multiplicities by mi = {#j : kj — i} — dimiP _1,1 (£). 
Thus S is given by m kl equations F klt i, . . . , F kl>mi of order ki, . . . , m kr 
equations F krA , . . . , F kr>mr of order k r . 

r 

We shall write 8 symbolically as Yl F ki = rriiEi, and call the latter 

i=i 

the type of £ (implicit: uj — 1 common characteristic and compatibil- 
ity). The following table shows types for systems of order /c max < 5. 



k 

| h 1 


1 


2 


3 


4 


5 


1 


E 1 




2 




2£ 2 


E 2 + -E3 

2£ 3 


£"2 + -£"4 
E 3 + £"4 
2£" 4 


£2 + E 5 
E 3 + £5 
£4 + E 5 
2E 5 












2E 3 + E 5 


3 






SE 3 


2£"3 + E4 
E 3 + 2£ 4 
3£" 4 


E 3 + E 4 + £5 
2i?4 + £5 
£3 + %E 5 

E4 + 2i?5 

3E 5 












3E4 + £5 


4 








AE A 


2£ 4 + 2E 5 

E4 + 3£/5 


5 










5£ 5 



The general (infinite) table is upper-triangular and the first row con- 
sists of one element only. 

The generalized Laplace transformations are arrows between the 
species in this table, which are directed toward the top-left corner. 
Here are sample routes for successive Laplace transformations: 



— y E2 + Eq — y E2 + E§ — y E2 + E4 — y E 2 + E 3 — y 2E2 — y E\ 
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and more complicated: 

> 2E 4 ->• 2E 4 + E 5 -> 3£ 4 -> 4£ 4 -> £ 3 + 2E 4 ->■ 

— >• 2i^3 4- £4 — >• 3 £-3 — ^ i?2 4" £3 — y E2 — y E\. 

We will discuss the arrows in more details in Section 12.41 

2.2. Complexity. The type notations introduced above suppress a 
lot of information about £ (branching into sub-cases), like its symbol 
class, the order of compatibility conditions etc. We shall introduce an 
important integ combination of these data, characterizing the 

complexity of a compatible system £ of class u = 1. 

Let Qi be the symbols of £ (see Appendix|B]for details). Starting from 
some jet-level the dimensions of these subspaces stabilize: dim^j = 1 
for % 3> 1. 

00 

Definition 1. Complexity of £ is the number x = ^(dim^j — 1). 
The following inequality holds for class u — 1 compatible systems: 

T h ^ ^min — ^max- 

r 

Proposition3. We have: x< {k^ — i) + {k\ — r) • (k r — r) . Equality 

i=i 

is attained for the 'boundary' cases: r = 2 and r = fc mm (=> r = k max ). 
Proof. We start with the case r — k\. Then letting k r+ i = 00 we get: 

dim - 1 = / *' i < h 

im9i \ (h-j), kj <i< k j+1 

and so the complexity equals 

kiCki — 1) , w, ,\ , r(r+l) 

i=i 1 

as required. 

For r < fci , x achieves the maximum value if the system Yli=i 
has one characteristic X only on the last step s = r, and has more 
before this; otherwise dim^j decreases more rapidly. In other words 
reducing the symbolic system by X, the system becomes of finite type 
only at the highest possible order. 

Thus the maximal growth of dim gi is achieved for 

{i, i < k\ 

(Ari - j), kj <i < k j+1 (j < r) 
(ki — r — j), i = k r + j (0 < j < k x - r) 

and then the direct check yields the desired value of x. 
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Alternatively dimension of symbols does not change if we assume 
that the system £ at order fc max = k r has more than 1 characteristic 
X, but we add equations of orders k' r+t = k r + t, 1 < t < k\ — r, so 
that only at the last step a unique characteristic is left. Thus we can 
use calculation with r = k\ since {k[ — i) = k r — r for r < i < k\. 

Finally if r = 2, then the reduced by X symbolic system is of com- 
plete intersection type and the equality follows from calculation of the 
dimension of the solutions space as in [KL^ . □ 



Thus complexity x is bounded via the orders of the system. On the 
other hand the amount of types grows super-polynomially with x. 

2.3. How many system types have the same complexity x? Let 

us list all systems of low complexity x = ^(dimgfj — 1) < 6. 

0. This is possible only for E\. 

1. Obviously the type 2E 2 . 

2. Only one decomposition (no intermediate zeros are possible) 
2 = 1 + 1. E 2 + E 3 . 

3. Two decompositions 3 = 1 + 1 + 1 = 1 + 2, we get resp. the 
types E 2 + E 4 and 3E 3 . 

4. Two decompositions 4 = 1 + 1 + 1 + 1 = 1 + 2 + 1, but the latter 
splits: E 2 + E 5 , 2E 3 and 2E 3 + £4 (two cases distinguish by no 
common characteristic and one common characteristic for 2E 3 ) . 

5. Three decompositions 5 = 1 + 1 + 1 + 1 + 1 = 1 + 2 + 1 + 1 = 
1 + 2 + 2. Corresponding types are: E 2 + E 6 , 2E 3 + E 5} E 3 + 2E 4 . 

6. Here we have four decompositions 6 = 1 + 1 + 1 + 1 + 1 + 1 = 
1 + 2 + 1 + 1 + 1 = 1 + 2 + 2+1 = 1 + 2 + 3, the third splits: 
E 2 + Ef, 2E 3 + Eq, E 3 + E4, E 3 + E4 + £"5, AE4. 

Denote by R(n) the number of u = 1 different types m iEi having 
complexity x = n. We have the following values for this function: 



x 
x 

X 
X 
X 

X 
X 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




asymptotics 


R(n) 


1 


1 


2 


3 


3 


5 


6 


9 


11 


13 




exp(ir\ / X n) 



One can prove (via a relation with the number-theoretic partition 
functions) that the quantity A defined by v^A = - ■ lim log R(n)/y/n 
satisfies | < A < |. 

2.4. Classification of transformations: Phenomenology. Let us 

summarize our results of §1.1111.41 Given a compatible linear system S 
of class oj = 1 we have constructed via generalized Laplace transforma- 
tion L : u i—)- Xu a new compatible system £ . We have observed that 
only three different situations were possible: 

(1) The transformed system £ is of class u = 1 and the inverse L _1 
is a differential operator. 
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(2) The transformed system £ is of class u = 1, but the inverse 
operator L^ 1 is obtained by solving a finite type system. 

(3) The transformed system £ is of class u = 0, but the inverse 
operator L _1 is integral. 

Note that situation (1) is generic. Finding the inverse in situation 
(2) is equivalent to solving a system of ODEs, while in situation (3) it 
is given by a parametrized ODE, so that the general solution depends 
on one unknown function. 

We claim that this is the general pattern, namely we have 

Theorem 4. For compatible linear systems of class u = 1 the general- 
ized Laplace transformations can branch into three situations (1), (2) 
and (3) above. The first occasion is generic. □ 

This will be demonstrated in the next section by showing that Laplace 
transformations decrease the complexity. 

We take by definition all systems of class u = to be of lower 
complexity than the systems of class u — 1. Thus situation (3), when 
we leave the table u = 1 satisfies the claim. 

Iteration of Laplace transformations leads either to a system of class 
lo = or to one equation Ei of class lo — 1. Both are reduced to the 



solution of ODEs in accordance with Sophus Lie's theorem [L[ . 

But for the linear u = 1 case and generic £ we obtained an algorithm 
for finding the solutions involving differentiations and quadratures only 
(for non-generic £ solutions of ODEs can occurs). 

3. Laplace transformations for linear systems 
Now after lots of examples we give a rigorous definition and proof. 

3.1. Generalized Laplace transformations. Consider a system £ = 
{F[u} = 0} of class lo = 1 on the unknown u = u(x,y) generated by 
a vector linear differential operator F = . . . , F r ). We can choose 
generators of £ in such a way that the maximal degree of Y in symb(Ffc) 
strictly increases with k (this is independent of basic/auxiliary gauges). 

Then there exists precisely one change X i— )■ X + a, a e C°°(M) 
such that in the decomposition Fk = ct k jjY^X l the number max{j : 
3kaj 7^ 0} is minimal (maximal j corresponds to the minimal k). 
This fixes the basic gauge, but leaves an auxiliary gauge freedom. 

Definition 2. With X fixed as above the generalized Laplace transfor- 
mation is the substitution 

u i — y v — Xu. 



Let 3(£) = [[Fx, . . . , F r }} be the differential ideal of £. 
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In order to define the transformed PDE system £ let us consider an 
ideafl J = {E i>0 Pij^X* G J(£)}, which factorizes J = J' ■ X for 
J' = {^2i>o ftijYiX 1 * 1 }. Let G = (Gi, . . . ,Gi) be a vector generator 
of J'. Then we define £ = {G[v} = 0}. 

Inverse Laplace transformation is defined as a scalar (but apriori not 
necessarily differential) operator L such that L ■ X — 1 mod 3(£). 

Proposition 5. Inverse operator is unique modulo f 3(£). 

Proof. Let L\ and L 2 be two inversions for the operator X, i.e. 

1 - L X X, 1 - L 2 X G =^ -L 2 )X e 3(f). 

This means that (Li — L 2 ) G or equivalently (Li — L 2 ) G □ 

Next let us consider the existence part. 

Theorem 6. For a generic system £ the generalized Laplace transfor- 
mation v = Xu has a differential inverse operator u = L[v]. 

Proof. Let us prolong £ to the place, where it becomes involutive. This 
is the jet-level k such that dimg% = 1. Discard all equations of order 
lower than k. Then we have r = k compatible PDEs of order k (type 
kE k ), so that x = ^ fc ~^ fc . 

The symbols of these equations are X k , YX k ~ l , Y 2 X k ~ 2 , . . . , Y k ~~ l X. 
Compatibility forces the first (k — 1) equations to be free of Y k ~ x terms. 
When the basic gauge is applied, the last PDE has the same property 
and then the compatibility implies that the first (k — 1) equations 
contain no Y k ~ 2 terms. 

This is symbolically shown in the following Young diagram (for 4E4), 
where we omit all terms in equations Fj of £ except Y l : 



> k 




fc-1 



It corresponds to the linear system of k equations and (k — 1) unknowns 
Y k ~ 2 u, . . . , Yu, u. Generically coefficients are such that the rank is 
maximal. Then we can exclude Y k ~ 2 u, . . . , Yu and get an expression 
u = L[v], where L is a differential operator of ord < k — 1. □ 

J is surely a left differential ideal, but if we restrict Y by [AT, Y] = 0, it will be 
also a right differential ideal. 
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More refined Young diagrams based on the type of symbol of £ can 
be drawn and then an inverse L with the minimal possible order (order 
of the inverse) can be chosen. 

For non-generic systems the Laplace transformation can have inverse, 
which is either a compatible Frobenius system or a parametrized ODE. 
Define order of the inverse in these cases to be or —1 respectively. 

3.2. Proof of the main result. Theorem [1] follows from 

Theorem 7. Under generalized Laplace transformation the complexity 
strictly decreases. Generically it decreases only by 1: x \- > x — 1. 

Proof. By re-scaling the unknown function u i— >■ a ■ u we can achieve 
X = d x , and we adopt this convention in the proof. 

Let us interpret complexity x via Cauchy data, namely a general 
solution of £ depends on one function of 1 variable, its q derivatives 
and on (x — q) constants. This follows from Cartan-Kahler theorem, 
and generically q = x. Moreover these are linear superpositions 

u = a q f {q) H h aif + a f + c x §\ H V c x - q (f) x - q , 

where aij = ai(x,y), <pj = <f)j(x,y) are some fixed functions, / = f(y) 
is an arbitrary function and are arbitrary constants. 

It is easy to check that a q depends only on y (this was the reason for 
the basic gauge as it is equivalent to u i-> u/a q ), so that the transfor- 
mation v — u x reduces x (only by 1 if {a. q -\) x ^ and (<f>i) x 7^ 0). □ 

Generically only transformations (1) from §2.41 are used and so the 
general solution is obtained from the composition of inverse Laplace 
operators u = A[f], where / is an arbitrary function of 1 argument 
and A a linear differential operator of order x. 

Transformations (2) with L solving Frobenius system of order 1 and 
transformations (3) from §2.41 with reduced system of order 1 and class 
u = (both cases belong to a generic stratum of the space of singular 
£) are equally good for Theorem [lj because linear finite type scalar 
systems of order 1 are solvable in quadratures. 

On the other hand, finite type systems of higher order are usually 
non-integrable in quadratures (and a theorem of S.Lie about solvable 
symmetry groups indicates when such integration is possible). 

For example, if the reduced system £ (or the system for inverse L) 
has type E\ + E 2 and class u = 0: 

f X 2 u + a 1 Xu + b l u = 0, 
1 Yu + a 2 Xu + b 2 u = 0, 

then it is generically non-solvable in quadratures. Indeed the first PDE 
is equivalent to a Riccati equation via a substitution u = exp J zdx. 
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Finally note that the algorithm for the Laplace transformation from 
Section I3TTI together with the first part of the proof of Theorem [7] imply 
that we can perform Laplace transformations so long the reduction has 
type u) — 1. Theorem H] follows. 

3.3. The role of relative invariants. The relative invariants found 
in the cases of §l.lH1.4l play the same fundamental role as the Laplace 
invariants, and they govern the branching into sub-cases of the types 
of £ , determining kind (1), (2) or (3) of the Laplace transformation. 

On the other hand we know that these kinds specify the form of the 
general solution, namely the amount (x—q) of constants it contains. 
These constants are the first integrals of the system £ and so can be 
detected by the internal geometry. 

Thus the generalized Laplace invariants of u — 1 type systems are 
the obstructions to existence of the first integrals. 

There can be one such differential invariant (as in the case T23) 
or several (as in the case T33), and there can also occur conditional 
invariants of various depth (that are defined only on the stratum given 
by vanishing of the previous invariants, as in the case T| 3 ). 

In our approach we first applied Laplace transformations to achieve 
decrease of the number q of the involved derivatives, and then had to 
find a (x — g)-parametric solution of u = type system. 

But equally well we could first restrict to the level leaf of the first 
integrals (freezing all the constants) and then solve a generic u = 1 
type system. 



4. Towards nonlinear theory 

4.1. Linearization theorem. From the internal viewpoint the system 
£ C J k (R 2 ), k = A; max , is a submanifold endowed with the induced 
distribution = C D T£, where C is the canonical Cartan distribution 
on the space of jets. 

We assume £ involutive, otherwise we need to prolong it to the jet 
level k such that dimgfc = 1. Then the distribution Cs has rank 3 - 
it is generated by the vector fields V x , V y (restricted total derivatives) 
and the vertical field (kernel of the projection T£ — > TJ k ~ 1 ). 

The first of these fields (or its combination with the third field) is a 
lift of our characteristic X and it coincides with the Cauchy characteris- 
tic X for this rank 3 distribution. Denote the (local) quotient manifold 
by £ = £/X and the quotient rank 2 distribution by LT = Cs/X. 

Let Iloo be the bracket- closure of the (generally non-holonomic) dis- 
tribution IT. We assume it is regular. The next claim is obvious. 
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Proposition 8. The amount of constants in the general solution of £ 
equals q = 00011111(1100), which also coincides with the codimension of 
the bracket- closure of Cs in £ . □ 

A regular distribution Ai = II generates the following: the weak de- 
rived flag Aj via the commutator of sections r(A i+1 ) = [T(Ai), T(Aj)] 
and the strong derived flag via r(Vj+i) = [T(Vj), r(Vj)], Vi = A x . 
The sequences of their ranks are called the weak and strong growth 
vectors respectively. 

Goursat distribution is the canonical rank 2 Cartan distribution of 
the jet-space J k (R) (we exclude singularities). If the distribution is not 
totally-nonholonomic, we can restrict to the leaves of the foliation of 
its bracket- closure. Distribution will be called Goursat-Frobenius if it 
is Goursat in all leaves, i.e. locally the distribution C x on J k ~R x M. m . 

Theorem 9. An involutive ui — 1 type system is internally (micro- 
locally) linearizable if and only if the corresponding rank 2 (regular) 
distribution II is Goursat-Frobenius. 

Proof. It is straightforward to check that both weak and strong flags 
grow in dimension by 1, so they are Goursat-Frobenius. The opposite 
direction is given by Cartan-von Weber theorem, which states that a 
rank 2 distribution has locally Goursat normal form if and only if both 
weak and strong growth vectors are (2, 3, 4, 5, ... ). □ 

As for external linearization, the responsible differential invariants 
are known only in some partial cases. The classical case of the 2nd 



order scalar ODE is due to S.Lie and R. Liouville [L2]. For higher 
order scalar ODEs the contact trivialization result is known |Dou] , and 
for 3rd and 4th order ODEs the linearization is done ((IMS, SMI] and 
the references therein), but to our knowledge no general linearization 
criterion is known (even for u> = type). 

For cu = 1 systems we expect our generalized Laplace invariants 
(which exist also in the non-linear situation) to play an important role 
in this classification problem. 

Corollary 10. If a PDE system £ ofu — 1 type is effectively lineariz- 
able, then £ possesses a closed form of the general solution. 

Here 'effectively' stays for algorithmic computability of the lineariza- 
tion transformation. This means the following. A Goursat distribu- 
tion can be transformed to the normal form by rectifying some vector 
fields. In the presence of symmetry algebra with nice (e.g. solvable) 
Lie structure this solution to ODEs can be made effective (for instance 
in quadratures). 
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4.2. Examples of nonlinear Laplace transformations. 

1. Consider the Liouville equation together with one of its Laplace 
integrals 

V'xx "2^x1 ^xy ^ ■ (6) 

This quasi-linear system is compatible of class u — 1, and it possesses 
no intermediate integral. Laplace transformation here is the same as 
in the linear case v = u x . 

The first PDE of gives the transformed equation E%: 

The second equation of yields the inversion: u = logv y . 

The solution v = — 2/(x + ip(y)) to E\ provides the solution to <Q 
via inversion: 

(x + ip{y)Y 

Notice that ([H]) is not linearizable by an external transformation of 
Moutard type (preserving x,y), but substitution w = —2m" 1 — x + 
2xe u u~ 2 makes an internal diffeomorphism of this system with the lin- 
ear 2 E2 

w xx = 0, w xy = (w y - w x )/x. (7) 
However according to E.Cartan |Ci| a Lie-Backlund type theorem holds 
in this case, so that every internal transformation is induced by a con- 
tact transformation. We can modify our external transformation to the 
following contact equivalence between §6§ and ((7]): 

u x 1 2e u u v 

x ^ — — , y H> y, w H> w = -x , u x i-)> w x = u y h> w y = —. 

2e u u x u% ■ u x 

In fact, one can show that ([7]) is the normal form for linearizable 2E 2 
without intermediate integrals (in the latter case we get u xx = u xy = 0). 

Remark. The standard Backlund transformation that linearizes the 
Liouville equation maps (T^j into non-linear system v xx = \v 2 x , v xy = 0. 

2. Consider another quasi-linear system 



u x _ y/U x U y 

VlXX " , ) V X y / 

x + y x + y 

for which the second equation is also Darboux integrable. 
Substitution u y = Q 2 transforms ([8]) into 

n = 2 ® x O = + - f9l 

x + y x + y [x + y) z 

This is a linear 2E 2 and it has Laplace transform v = Q x to 

2v 
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with the inverse Q = (x + y) 2 v y + (x + y)v. But the lift (jHJ) to (jEJ) is 
given by the (compatible) Frobenius system 

u x = (x + y) 2 Q 2 x , u y = Q 2 . 

Another nonlinear Laplace transform is given by the following dif- 
ferential substitution with differential inverse: 

w 2 

u = w — xw x — ^ w = u + (x + y)u x + x(x + y)\fuZ- 

x + y 

This transforms (JSJ) into another quasi-linear system 2E 2 

W XX = 0, Wly + XW X y ~Wy = 0. 

This readily yields the general solution u = My) — <p'(y) = ip'(y) 2 

x-\-y 

of (jHJ), and consequently its closed (non-Moutard) form solution: 

2 a t \ 1/ \ t \ (tcr"(t) - a'(t)) 2 
x = s, y = a"{t), u = tV'(t) - 2ta\t) + 2a(t) - 1 ^ — f )} ■ 

3. Consider fully non-linear system of type 2E2 

3U XX V,l y = -l, U X y=-l/Uyy, (l0) 

where the first equation belongs to the Goursat's examples of Darboux 
integrable systems [G] . 

This case is more complicated and the nonlinear Laplace transfor- 
mation is given by A = u 2 y (notice it is of the 2nd order). The result 
of the transformation has type E\ and is the equation of gas dynamics 

\ x = \ Xy. 

The inverse is however not differential, but is given by the Frobenius 
system (compatible 3E 2 of class u = 0): 

U> XX gA ^ , ^ X y A ^ , ^yy ^ ^ ; 

and thus (TlOT) is solvable by quadratures. Notice that the choice of A 
corresponds to the characteristic of (TTUI) . 

4.3. Examples of nonlinear Laplace invariants. Consider 

Vi ^> H'xi U"yj "in ^ X yi V"yy) 0; 
G\Xi y, U, U x , Uy, U XX ] U X y, Uyy) 0, 

Linearization of this system has the form 

X 2 U + a x XU + b x YU + c x U = 0, 
YXU + a 2 XU + b 2 YU + c 2 U = 0, 

Compatibility implies b\ = 0. Notice that this does not follow directly 
from the arguments for the linear case, because the bracket of lineariza- 



£ 



tions differs from the linearization of the bracket |Ki| . However we can 
use the following straightforward statement: 
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Proposition 11. For nonlinear operator F e diff and an operator L 
in total derivatives we have: £l-f — L ■ £ F = mod (V a F). □ 

Now the system {F = G = 0} for the operators of the above form 
is compatible iff there does not exist operators in total derivatives 
(C-differential operators) Li,L 2 such that £l 1 -f+l 2 -g £ (1 3 X, Y, Y 2 ), 
which, by the above proposition, implies L\ ■ tp + L 2 • £g G (1, X, Y, Y 2 ) 
on the equation £(°°). And this would imply b\ — and in the case 
b 2 = also c\ = for the linearized system 

The gauge X = X + b 2 transforms the system to 

J X^ + cuXtZ + Cif/ = 0, 
\ FXf/ + a 2 XU + c 2 U = 0, 

where a\ = a\ — 2&i, ci = Ci — a^b 2 — X{b 2 ) = (due to compatibility), 
a 2 = a 2 and c 2 = c 2 — a 2 b 2 — Y{b 2 ) (notice similarity of the expressions 
for c\,c 2 with the classical Laplace invariants, see Appendix [A]) . 

This function c 2 is a relative invariant with respect to gauge transfor- 
mations. It will be called nonlinear Laplace invariant of the problem. 
It vanishes precisely when there exists a (higher) intermediate inte- 
gral / for the system, i.e a function I on jets with the vanishing total 
differential d(I) = due to the system £. 

In our case c 2 = implies existence of an integral J. Indeed, the 
substitution V = XU reduces is to the Frobenius system 

(X + a 1 )V = 0, (Y + a 2 )V = 

with an intermediate integral <&(x,y,V) = c, which can be integrated 
to give the integral <j)(x, y, Vu) = c for the original system. 

Example 1. The Liouville equation coupled with I^-intermediate in- 
tegral is a compatible 2E 2 of class u — 1: 

^xx 2 ^x t^xy ^ 0* 

Its linearization writes 

X 2 U - u x XU = 0, YXU - e u U = 0. 

Already in this form b\ = b 2 = c\ = and so c 2 = e u ^ 0. Conse- 
quently this pair (and the Liouville equation alone) has no intermediate 
integrals. 

Example 2. Now consider such 2E 2 : 

^xx U>x& •} ^xy UyC . 

with linearization 



X 2 U - e u XU - u x e u U = 0, YXU - e u YU - u y e u U = 0. 
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In accordance with the above theory b\ = 0, but b 2 = —e u ^ (starting 
from these coefficients we can calculate by the above formula c 2 = 0). 
So we make a change of frame X = X — e u and get the system 

XXV = 0, YXU = 0, 

with an obvious intermediate integral XU = c. This linearization can 
be integrated to obtain the first integral of the original 2E 2 : I — u x —e u . 

Now let us turn to the case E 2 + E 3 , for which the theory is similar. 
Example 3. Consider the following linear ui — 1 system of type 

U xx l^xyy •E'U'xy ~t~ U y 0* 

Laplace transformation v — u x maps this to E\ : v x — 0, and the 
inverse can be found from the Frobenius system 

UX = V, Uy = XVy - Vyy. (ll) 

Integration yields the intermediate integral: / = u xy — xu x + u = c. 

Remark. Existence of integral I here does not follow from vanishing of 
a relative differential invariant: we have 4 / in the notations ofT^. 
But the symbol of this system is degenerate, which means vanishing of 
certain relative algebraic invariant. 

Example 4. Now let us consider a nonlinear example: 

1 3 

u X x g^j/y ^ x yy ^yy^yyy 

Its linearization is 

V 2 X U - u 2 yy V 2 y U = 0, V x V 2 y U - u y yVp - u yyy V 2 y U = 0. 

Let us introduce the operators X = T> x — u yy T> y (characteristic) and 
Y — T> x + UyyVy. They do not commute [X, Y] = u yyy (Y — X). Notice 
also that X{u yy ) = and Y{u yy ) = 2u yy u y y y and using this we re- write 
the above linearized E 2 + E 3 system so: 

XYU = 0, X 3 U = 0. 

Notice that b 2 = d 2 = e 2 = in the notation of T23, but d\ — only 
in this form and d\ 7^ in another representation of the first equation 
of the system: 

YXU + Uyyy(YU - X U) = 0. 

But the coefficient e± in this normal form is a relative differential invari- 
ant and it controls existence of the integral. For the linearized equation 
it equals 

-^X 2 U = c, 

Uyy 
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which can be re-written (modulo equation) as V x V y U — u yy VyU = c, 
and this can be integrated to obtain the intermediate integral of the 
original problem: 

I — u xy — 2 U yy 

We get the following u — 1 system 2E 2 that is compatible for any c 

13 12 
u xx — 3 U yyi U xy = 2 U yy C ' 

5. Conclusion: Beyond u — 1 and other generalizations 

The theory of Laplace transformations developed for linear 2nd order 
PDEs (u = 2) was further modified to work in the non-linear situation 
by Darboux and Goursat. 

In this paper we mainly concentrated on the general linear u = 1 
theory, leaving the non-linear case for a separate publication. However 
even the linear theory sheds a light on the obstructions for u> > 2. And 
indeed this latter case is poor compared to the theory of u = 2 type. 

An important effort to generalize the classical (linear) Laplace theory 
to operators of higher order was undertaken in [R]. In particular it was 
observed that the Laplace transformation rule applied to one equation 
inevitably generates several equations, which was the main obstruction 
for effective integration theory. 

Another attempt was jP], where transformation theory was con- 
structed for some particular class of equations, with a degenerate sym- 
bol. The reason for this is that starting from 3rd order the scalar 
PDEs on the plane have invariants of the symbol, since the latter can 
be considered as a planar web. Generic webs have lots of independent 
differential invariants, and they admit no symmetries. 

For recent advances of the factorization part we refer to [Tj and to 
vast range of papers on differential Galois theory. But the transforma- 
tion part has not seen much progress beside the classical 2nd order case 
(we refer for the development of the Darboux theory to |AF|, IA IAK| ) . 
The reason is that the Laplace transformation for u> > 2 inevitably 
increases the complexity x and this gives less chances of termination 
for the sequence of Laplace invariants. 

Laplace theory, extended for systems of second order equations as the 
theory of multidimensional conjugate nets in jD], was further developed 
in |KTl IFej . This worked for general u but again for equations of very 
special type, namely semi-Hamiltonian and integrable. 

Thus the general theory for u > 3 is still lacking, and below we 
briefly discuss these cases and u — 2. The case u — 1 turns out, on 
the contrary, to be a perfect arena for Laplace ideas. 
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5.1. Generalized Laplace transformations for uj = 2. The classi- 
cal method of Laplace concerns E 2 , we recall it in Appendix |A] The 
important distinction of this case from u = 1 is that we have a pair of 
Laplace transformations v = Xu and w = Yu. 

Consider the case 2E 3 . A system £ of this type with generic symbol 
writes 

j YX 2 u + ai X 2 u + hY Xu + Cl Y 2 u + d x Xu + e x Yu + f x u = 0, 
\ Y 2 Xu + a 2 X 2 u + b 2 YXu + c 2 Y 2 u + d 2 Xu + e 2 Yu + f 2 u = 0. 

Compatibility implies c\ = 0. We begin with Laplace transformation 
via characteristic X. Then the basic gauge yields c 2 = and compati- 
bility implies e x — 0. 

With these changes we let v = Xu to be the direct transformation 
and find the inverse u = L\V from the 1st PDE of £, where Li = 
—/^(YX + a\X + b{Y + di) and we assume f\ ^ 0. 

Inserting the inversion to the substitution gives the first transformed 
equation {XL\ — l)v = 0, and the second is obtained from the 2nd 
PDE of £: (e 2 Y + f 2 )u = L 2 v, where L 2 = —(Y 2 + a 2 X + b 2 Y + d 2 ), 
namely {{e 2 Y + f 2 )Li — L 2 )v = 0. The transformed system has symbols 
YX 2 ,Y 2 X and the type 2E 3 unless e 2 = 0. In the latter case the 
transformed equation is of type E 2 . 

When fx — 0, the compatibility gives e 2 = 0. Provided / 2 ^ 
we get inversion from the 2nd PDE of £, and the 1st PDE gives the 
transformed £ = {(YX + a x X + b x Y + d x )v = 0} of type E 2 . 

In the case f\ = f 2 = the transformation v = Xu brings £ to a 
system 2E 2 of type u = 1 and the inverse m = X _1 f is integral. 

Situation with the second transformation w = Yx is similar. 

Definition 3. Complexity of uj = 2 type system is defined by the for- 
mula 



The classical case E 2 has x = and 2E 3 corresponds to x = 1 
We can also study the cases of higher complexity E% + E4 (x = 2), 
E 3 + E 5 , 3E4 (both x = 3) etc. forming the zoo of uj = 2 type. These 
cases are similar and following the considered pattern we get such a 
conclusion. 

Theorem 12. The generalized Laplace transform (any of two) foruj = 
2 type generically preserves the type and complexity. In singular cases 
it can decrease x ; but leave uj or it can reduce the class to uj = 1. □ 

Note that we cannot reach the systems of finite type uj = in one 
step as our Laplace transformations are always of order 1 (higher order 
generalizations are plausible). 
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The integration theory in the general case co = 2 is thus characterized 
by dropping to class u> — 1 type in a sequence of Laplace transforma- 
tions to both ends. This would reduce solution of £ to ODEs. 

If only one sequence of transformations (say v = Xu) terminates at 
ou — 1 class, this gives semi-integrability, i.e. possibility to find many 
solutions (a family depending on 1 function of 1 variable) but not all. 

5.2. Generalized Laplace transformations for cu > 3. For the 

class cu > 3 we define the complexity by the formula 

oo oo 

x = y^(dim gi — cu) + §tu(cu — 1) = ^^(dim g-i — cu). 

i=0 i=ui 

The simplest cu = 3 case, corresponding to x = 0, is of type E 3 : 

ZYXu + aX 2 u + bYXu + cY 2 u + dXu + eYu + fu = 0. 

Here Z is a linear combination of X and Y. If the web of characteristics 
is parallelizable, we can assume Z = X+Y. Now there can be 3 Laplace 
transformations, with symbol X, Y and Z respectively. 

Let us study the X-transformation. By the basic gauge we achieve 
c = 0. We assume for simplicity that (the new coefficient) e ^ 0. 

Then the Laplace transformation v = Xu leads to = Lv, where 
^ = Y + £,L = ZY + ... This yields the equation (XL -V)v = u, 
where 6 = [X, \&] is a function. 

Only in the singular case 6 = has the transformed equation type 
E3 again: then £ is (XL — ty)v = with the same symbol ZYX. 

In the general case 9 7^ the above formula gives the differential 
inversion u = 9^ 1 (XL — ty)v and inserting this into the equation and 
the substitution we get the following pair of PDEs: 

(ZYX 2 + . . . )v = 0, (ZY 2 X + ..> = 0. 

Thus the result of the transformation is a system £ of type 2E^ (and 
class 00 = 3). After the next transformation we get equations with 
higher complexity x, like 3E 5 or E4 + E 5 etc. 

Theorem 13. For to = 3 type system any (of three) generalized Laplace 
transformation generically increases the complexity by 1. Only in some 
singular cases it preserves the type or decreases x or reduces 00. □ 

For general u the complexity increases by < u — 2 units with any 
generalized Laplace transformation. Consequently the integration in 
closed form is a rare occasion for the class 00 > 3. 
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Appendix A. The method of Laplace 

This method is well-known [D|, [Gj |F] . But for completeness, and to 
show parallel with what is done for u — 1 case, we give a short review. 

The classical case concerns a hyperbolic equation £ of the type E 2 
(which has class oj — 2). Consider at first the linear case: 

A[u] = u xy + au x + bu y + cu = 0, 

where a, b, c are functions on M = R 2 . We try to factorize 

A[u] = (V y + a){u x + bu) — k u 

with k = by + ab — c the 0th level Laplace invariant^. If k = 0, we 
have an intermediate integral v = u x + bu. 

Otherwise the latter formula defines the Laplace transformation with 
differential inverse u = -j^(v y + av). Substitution of the inverse to the 
direct Laplace transformation gives the equation 

A x [v] = (D x + b) (^^) - v = v xy + a x v x + b x v y + c x v = 

Since the equations are effectively equivalent, we search for factoriza- 
tion of this new PDE. The obstruction for existence of intermediate 
integral is another relative invariant k\ = {b\) y + a\b\ — c±. 

If the sequence k ,ki, . . . stops at zero k n = 0, £ is called Darboux 
semi-integrable. In this case the equation £ possesses an intermedi- 
ate integral of order n, which together with the original PDE forms a 
compatible pair E 2 + E n . Thus semi-integrability can be interpreted 
as a reduction of = 2 equation to a system of class 00 — 1, which is 
integrable via ODEs. 

The method of Laplace assumes that the sequence of invariants is 
finite in both sides, where to the other side we add up invariant^ 
ho, hi, . . . obtained via X-factorization 

A[u] = iV x + b){u x + au) - h u, 

so that the corresponding Laplace transformation is v = u x + au with 
differential inverse provided h = a x + ab — c 7^ etc. 

Remark. The sequence . . . hi, ho, k , ki . . . of relative invariants leads 
to a collection of absolute invariants, but they do not form a basis that 
solves the equivalence problem; some other invariants shall be added [I]. 



It belongs to a sequence of (fundamental) relative invariants with respect to 
linear point transformations preserving the symbol of £ . 

9 The invariant h n can be also interpreted as an obstruction to find a compatible 
differential constraint for the hyperbolic equation £ = {A[u] = 0}, which is a 
y-parametrized ODE of order n: J2o otiU^ 1 ^ = 0. 
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In the non-linear case we cannot proceed with precise transforma- 
tions, but can compute the sequence of Laplace invariants anyway via 
linearizations of the corresponding operators \AK\ IAJ] . 

Darboux integrability, i.e. vanishing of these generalized Laplace in- 
variants to both sides, again implies closed form of the general solution, 
as in the linear case. The above interpretation of semi-integrability via 
reduction of the class from u = 2 to u = 1 is still valid. 



Appendix B. Spencer 5-complex and cohomology 

These will be described only in the linear case and for the base 
M = M. 2 , see the general description in |KLi| . 

In the case of scalar equations the symbols g k are the subspaces of 
the kernels 7ik,k-i '■ J k M — > J fc_1 M, which can be identified with S k T*, 
T = T X M being the (model) tangent space. 

The equation of the lower order S kl specifies the subbundle £ kl C 
J kl M with the fiber g kl . The next symbols g k +i equals the prolonga- 
tions = (g k <g> T*) n (S k+1 T*) for h < k < k 2 . 

They correspond to the symbols of the prolonged equation £ kl pro- 
vided the compatibility conditions of order < fc+1 hold. At the jet-level 
k 2 new equations are added, and we get g k2 C g k ^ C S k2 T* etc. 

These symbols are united into the Spencer 5-complex 

-> g k +i -> g k ® T* A g k _! ® A 2 T*, (12) 

where the first arrow is the inclusion and the morphism 5 is the Spencer 
differential, i.e. the symbolic exterior derivative (if we interpret g k ®T* 
as differential 1-forms on T with polynomial coefficients). 

The cohomology group H k,l (S) at the mid-term of f[T2"j) counts (i.e. 
dim if equals) the amount of new equations in £ at the jet-level k. 

The cohomology H k ~ 1,2 (S) at the last term of ([12]) is the most im- 
portant — its elements are the compatibility conditions. 

If all the compatibility vanish, the system £ is said to be formally 
integrable. This is tantamount to the claim that the projections ir k +i ,k '■ 
£ k -\-i — >■ S k are vector bundle^ (regularity means constancy of ranks). 

Formally integrable systems S are not necessarily involutive on the 
level k = fc max , but they are such after some prolongations (for larger k). 
Since involutivity is equivalent to vanishing of the Spencer cohomology 
jgj jjk,% _ h_ a pp ens on t ne jet-level k where dim g k stabilizes. 

The growth of dim.g k for large k is given by the Hilbert polynomial 
Ps(k) = cr k d + . . . , where d is dimension of the characteristic variety 



Char t '(^) C P C T* and a is its degree [KL; . 



'Contrary to this, solvability only means that £oo — > M is a (nontrivial) bundle. 
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For a scalar system £ the (complex) characteristic variety Char (£) 
is the intersection of all characteristic varieties for individual equations 
Fi = from £ = {Fi = : 1 < i < r} (for each of them the character- 
istics are defined as loci of complexified characteristic polynomial |Ptj . 
namely of the Fourier transform of the symbol smbl(Fj)). Alternatively 
a covector p G C T* \ {0} is characteristic if p k G for all k < k max . 

Notice that for linear systems the variety Char (£) depends only on 
the point x of the base M = M 2 . 

If £ is determined or overdetermined, then Char {£) consists of oj = 
a points corresponding to (complex) characteristics of £, i.e. d = 
and we obtain that = uj for large k. 
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